Charge Ordering Due to Magnetic Symmetry Breaking 
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It is argued that both transitions observed in 50%-doped manganites, at the Neel temperature 
(Tn) and the so-called charge ordering temperature (Tco), are magnetic. Tn corresponds to the 
order-disorder transition, which takes place between ferromagnetic zigzag chains, while the coherent 
motion of spins within the chains is destroyed only around Tco- The magnetic structure realized 
below Tco is highly anisotropic. It is dressed by the lattice distortion and leads to huge anisotropy 
of the electronic structure, which explains stability of this state as well as the form of charge-orbital 
pattern above Tn. The type of phase transition at Tn depends on lattice interactions. 

PACS numbers: 75.25+z, 75.30.Vn, 71.23.-k, 75.10.Lp 
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The so-called "charge-ordered" state of the colossal 
magnetoresistive perovskite manganites is one of the 
most interesting, important, and at the same time - 
somewhat obscured phenomena Q. Since the charge 
ordering (CO) is accompanied by the orbital ordering 
(00), and sometimes - by the very peculiar zigzag 
antiferromagnetic (AFM) ordering (Fig. [IJ, it always 
causes a problem with the identification of the true 
order parameter responsible for the physics of this 
complicated phase. Despite the long history, a visi- 
ble progress in this direction has been made only re- 
cently 0, 0, 0, H, @- It appeared that many low- 
temperature (T) properties of CO manganites, like 
Ndi^Ca^MnOa (NCMO), Ndi^Sr^MnOs (NSMO), 
and Lai_ K Sri +2: Mn04 (LSMO) at x=0.5, can be under- 
stood on the basis of purely magnetic considerations, by 
combining traditional ideas of the double exchange (DE) 
with the anisotropy of transfer interactions between 
two e g orbitals |8j, and formally without involving addi- 
tional mechanisms of the CO and OO. 



It is true that in order to stabilize the zigzag AFM 
state, one should first break the cubic symmetry. Most 
of attempts were devoted to the search of some exter- 
nal (with respect to the DE) mechanism, which would 
first stabilize the CO/OO pattern shown in Fig. and 
thereby predetermine the form of the zigzag AFM state. 
The main idea proposed in Refs. 0, LH was that in the 
degenerate DE model, the cubic symmetry can be bro- 
ken by internal means, by imposing the zigzag AFM or- 
der which forms the desired CO/OO ordering pattern, 
and this symmetry breaking is already sufficient to sta- 
bilize the zigzag AFM state. Thus, the appearance of 
the CO phase at T—0 can be explained by the following 
formula: anisotropic AFM order — > symmetry breaking 
— > stability of this AFM order. The last step is consid- 
erably facilitated by the fact that the zigzag AFM phase 
is a band insulator 0, 0, Q, ■ The finding is very im- 
portant, because the magnetic origin of the CO phase 
naturally explains the fact that it can be melted by the 




FIG. 1: Magnetic structures in the CO phase. Two AFM sub- 
lattices below Tn are shown by different colors. Above Tn, 
the Neel state is replaced by the inter-chain disorder, which 
coexists with nearly perfect FM order within the chains. New 
magnetic unit cell is denoted by x-y. The anisotropic mag- 
netic structures give rise to the OO patter, which is coupled 
with the oxygen displacements, shown by open circles. The 
CO is only illusive and reflected in the existence of two non- 
equivalent sites (A and B), produced by the magnetic sym- 
metry breaking. 



magnetic field, giving rise to negative colossal magnetore- 
sistance (CMR) effect The second mechanism, which 
additionally stabilizes the zigzag AFM state 0, Q , is the 
Jahn-Teller (JT) effect pjj. 

Is this the end of the story? Apparently not, because 
this picture is believed to be at odds with the tempera- 
ture behavior of CO manganites, which can be of two 
types. The first one, realized in NCMO and LSMO, 
is characterized by the distinct transition temperature, 



Tco, which is manifested in appearance of the CO/OO 
superlattice spots in the neutron [llj and x-ray [12j ex- 
periments, as well as the sharp increase of the resistiv- 
ity, and the peak of magnetic susceptibility Q, . The 
Neel temperature, Tn, is mainly manifested in the ap- 
pearance of AFM superlattice spots 0, . The typical 
values of T co and Tn, observed in LSMO, are 217 and 
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110 K, respectively. Since Tco>7n, the CO is typically 
regarded as a more fundamental effect, which drives the 
AFM transition at In. Another possibility, realized in 
NSMO, is the first-order transition to the metallic fer- 
romagnetic (FM) phase above Tn, which coincides with 
Tco- More generally, £=0.5 is the point of multi-critical 
behavior of CMR manganites |l3j ]. 

The purpose of this Letter is to argue that such a tem- 
perature behavior can be understood in terms of the par- 
tial spin disorder. Namely, only the long-range AFM 
order in the ^-direction disappears at Tn, whereas the 
FM coupling within the chains survives this transition. 
Thus, in the interval T^<T<Tqo we deal with the one- 
dimensional spin disorder in the direction perpendicu- 
lar to the chains. The disorder is highly anisotropic. It 
breaks the cubic symmetry, stabilizes the FM chains and 
gives rise to the CO/OO pattern even above Tn fnHlijl. 
We will concentrate on the dynamics of two-dimensional 
pattern shown in Fig. ^ 

The behavior of CO manganites at T=0 is well de- 
scribed by the following model, combining DE physics 
with the harmonic JT effect lal: 



TLij = 



(1) 
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[14j. e; is the spin direction at 
is the matrix of transfer interac- 



where ^ij = ^ 1+ ^' ej 

the site i. Uj = \\tij L „ 
tions between nearest- neighbor (NN) e g orbitals: x 2 -y 2 
(L,L'=1) and 3z 2 -r 2 (L,L'=2), in the local frame spec- 
ified by the oxygen sites (corresponding to the xy and 
3z 2 -r 2 orbitals in the xy-ir&me in Fig.^l. The NN hop- 
pings along the (±1,1) directions are parameterized ac- 
cording to Slater and Koster: tij = j(2^fy/3a z +a z ), where 
a x and a z are the Pauli matrices. Ai=||A^ L || describes 
the e g levels splitting caused by the cooperative JT 
distortion Sj T =(d L -d s )/{d L +d s ): A A =-^(V3a x +a z ), 
Aa>=^(V3&x— &z), and Ab——\o z (A being the 
electron-phonon coupling constant). The considered dis- 
tortion is compatible with the magnetic zigzag structure 
(Fig. 2]). £(<5jt) and A are estimated from Harrison's law 
for the p-d transfer interactions Q. This yields A~14to- 
J <0 is the AFM superexchange (SE) coupling between 
NN spins (\J s \~0.lt o H3) and 2Ci=£„ eN N e ; • e n . The 
spring constant k may vary from several tens to hundreds 
eV [lOj (below we will give a more realistic estimate for k 
from the analysis of the phase transition at Tn). to=t(0) 
is used as the energy unit (t =0-7 eV 0). 

What goes first in the Hamiltonian JTJ and causes the 
required anisotropy of the electronic structure: the mag- 
netism or the JT effect? In the harmonic approxima- 
tion around a homogeneous magnetic state, the ground 
state may have only one minimum at <5jt=0 jlfij . There- 
fore, the JT distortion alone cannot be the source of the 
anisotropy. On the other hand, the magnetic part of 



the Hamiltonian JQ) can, and the equilibrium magnetic 
structure at T=0 can be successfully described in terms of 
competition between DE and SE interactions 0, 0, 0, • 
Therefore, even for the finite-temperature description, it 
seems to be logical to start with purely magnetic Hamil- 
tonian by enforcing <5jt=0 in Eq. ft}. What do we expect 
in this case? 

The stability condition of the zigzag AFM phase at 
T=0 is given by J 1 AB < \ J s \ < jj^ , where J 1 J B (J A l B ) are 
the NN DE interactions within the same (between differ- 
ent) zigzag chains 0. For Sjt—0 the exchange anisotropy 
AJ=J AB — J A ^ B is of the order of 0.07 (and becomes 
larger for (5jt>0). Our basic idea is to connect this 
anisotropy, which roughly corresponds to the tempera- 
ture range of 550 K, with the existence of two transition 
points, Tn and Tco- Since | J s \ is close to J AB (about 
0.106 for <5jt=0), the AFM coupling between the chains, 
J^ AB +J S , will be much weaker than the FM interaction 
within the chain, J AB +J S . Therefore, the long-range 
AFM order in the x-direction will be destroyed first with 
the increase of T. This corresponds to the first transi- 
tion point, Tn, and yields the following model for the 
spin disorder below Tn- 

All spins in an arbitrary chain I are rigidly coupled 
by J AB (the effect of inter-chain disorder on J AB will 
be discussed below). This imposes the constraint on the 
spin directions within the chains: ei e g=e£. The distri- 
bution function for ei=(ef, ef, ef), which is required in 
order to calculate the thermal (or orientational) aver- 
ages, can be taken in the form P(e^)~exp(— hef), cor- 
responding to the mean-field approximation for inter- 
chain interactions. Then, £ij = l for i and j belonging 
to the same chain. Our goal is to calculate the elec- 
tronic Green function Q (corresponding to the first two 
terms in Hy), averaged over {e^}. As the zeroth order 
approximation along this line, all £y connecting neigh- 
boring zigzag chains are replaced by their orientational 
averages, £«±i— [jj. Corrections to this approxima- 
tions can be constructed in a systematic way by using 
the perturbation theory expansion with respect to fluctu- 
ations A£=£«-ti-— £ with subsequent averaging over {e^}. 

Then, Q has the following form, after transformation to 
the reciprocal space (k y ) along y and using matrix repre- 
sentation for the orbital, chain, and inequivalent atomic 

indices: % K (e)={l + £~ 2 Wr[G° ky y ] n } K ^ 
where Q% corresponds to the zeroth order DE Hamilto- 
nian, and ~5u±i is the super-matrix of transfer inter- 
actions between NN chains. In these (shorthanded) no- 
tations, A£'s have the same chain indices as tj; 's. Since 

A£=0, the first-order correction vanishes, and the matri- 
ces t i , which enters the matrix multiplication for higher 
n, should have at least one common chain index in order 

to contribute to Q. Then, it is feasible to go up to n—2 



3 



and 3. In most cases, n=2 is already sufficient ^t|- Fi- 
nally, Gk y (e) can be transformed to the real space <?ij(e). 

The free energy T is estimated using variational mean- 
field approach [7|. First, we calculate the DE energy: 



E 



DE 



de In 



1 + exp 



k B T 



Tr^e). 



with i running over the unit cell shown in Fig.^ £p be- 
ing the chemical potential, and Tr^ denoting the trace 
over L. Then, -Ede is combined with the SE energy 
and the entropy term, corresponding to the molecular 
field h 0- This yields T(h), which should be mini- 
mized with respect to h or the AFM order parameter 
m J 4=tanh _1 h—h^ 1 . 

The orbital polarization at the site A (the difference of 
atomic populations of the 3x 2 -r 2 and y 2 -z 2 orbitals, in 
the oxygen frame) is given by 

A = -^Jm J def T (e - e F )Tr L {(V3d x +d z )V A A(e)}, 

where /t(e _ £f) is the Fermi distribution function. (Sim- 
ilar parameter Ob related with As at the site B is of the 
order of 0.1, and only weakly depends on the degrees of 
the JT distortion and the spin disorder). The charge 
disproportionation 

C = -~Im J def T (s - e F )Tr L {d AA (e) - d BB (e)} 

shows the degree of CO. The stiffness of the FM zigzag 
chain is controlled by [lij 

J AB = ~^~ Im J d£ h{e ~ e F )TT L {t AB g B A(e)}. 

In order to be stable against the thermal fluctuations 
J^ B should satisfy the inequality: (J^ B +J s )^>ksT. 

Results for 5jt=0 are summarized in Fig. The or- 
dered AFM state is a band insulator HHEI3. Wil1 it; 
be insulating in the case of the spin disorder between the 
chains? In the zeroth order of A£, it is equivalent to the 
change of £«±i in the DE Hamiltonian from £=0 (m A —l) 
to £=| (m A =0) 7J. The corresponding electronic struc- 
ture is semi-metallic (SM), meaning that there is a di- 
rect gap in each k-point of the Brillouin zone. How- 
ever, a finite overlap between different k-points closes the 
real gap. Thus, the DE mechanism alone is insufficient 
to make the disordered zigzag phase a band insulator, 
though the situation is certainly close to it, and as we 
will see below the insulating behavior can be easily pro- 
duced by the JT distortion. 

The SM character of the electronic structure prede- 
termines the type of the phase transition at Tn. The 
collapse of the band gap at small m A leads to the gain of 
the DE energy and stabilizes an additional minimum of 
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FIG. 2: (a) Effects of inter-chain disorder on the electronic 
structure for <5jt=0. The disorder closes the band gap and 
leads to the semi-metallic (SM) behavior for small m A . (b) 
Free energy T against m A for T~Tn and |J S |=0.13. The 
SM electronic structure leads to the additional gain of the 
DE energy and stabilizes the new minimum at m A =0. (c) 
and (d) show the orbital polarization Oa against m A and T. 
The form of T{m A ) predetermines the first-order transition 
and discontinuity of O a (d). Dot-dashed lines, solid lines and 
symbols £g> correspond to the zeroth, second and third order 
of A£, respectively. Tn, obtained from T{Q)=J-{mA), is 0.15 
and 0.10 in the zeroth and second order of A£, respectively. 



T at m A —0, which coexists with another minimum cor- 
responding to an insulating phase at finite m A . Thus, the 
phase transition will be of the first order, which is charac- 
terized by the abrupt change of O a . C vanishes without 
distortion Q. This scenario is qualitatively consistent 
with the behavior of NSMO compounds. It also explains 
the observed phase coexistence around Tn (l8lll9| . which 
may lead to the percolation and the CMR effect [2fl |. 
However, most of CO manganites, like NCMO, exhibit a 
continuous second-order transition at Tn- 

What is missing in the above picture? Below, we will 
argue that the difference between NSMO and NCMO is 
controlled by the ratio A//s @, and can be understood by 
considering the response of the crystal structure to the 
highly anisotropic magnetic structure. First, we estimate 
the force /=— dEvE/dSjTls T=0 acting on the lattice de- 
grees of freedom and coming from the electron-phonon 
term in Eq. QJ 21], which couples (5jt with the parame- 
ters of orbital polarization. This yields (per one unit cell) 
f=-2X(0 A +0 B ), which is of the order of 10 eV. Then, 
we note that even a modest distortion 5jt^0.04 ji^] is 
sufficient to open the band gap in the disordered zigzag 
phase (Fig.|3J). Using /, we estimate the spring constant 
k which would correspond to such equilibrium distortion: 
«=— f /4J.TT- This yields k~65 eV, in agreement previous 
estimates |lCj . The temperature behavior corresponding 
to this ratio \/n is shown in Fig. Et- Since the system is 
insulating for all m A , J-(m A ) has only one minimum for 
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each T. Therefore, the transition at Tn will be of the sec- 
ond order. Finite <5jt enhanced the orbital polarization 
Oa and gives rise to the small charge disproportionation 
C. Both parameters are affected by the spin disorder 
and monotonously decrease with T, that explains sim- 
ilar intensity change of the CO/OO superlattice peaks 
in the neutron scattering experiment Tn, obtained 
from d 2 ^/dm\\ m Q —0, depends on the JT distortion. 
For (5j T =0.04 and |J S |=0.09, it is of the order of 120 K, 
which is close to the experimental value 110 K observed 
in LSMO QEl. Then, the ration is ( J^ B + J s )/k B T N 
is about 4, meaning that the FM coupling within the 
chains remains robust even around Tn. This justifies the 
proposed model of inter-chain disorder. 




FIG. 3: (a) The electronic structure of the disordered zigzag 
phase with and without JT distortion. Sjt opens the band 
gap. (b) Phase diagram showing the regions of semi-metallic 
(SM) and insulating (I) behavior versus <5jt and iua- (c) 
An example of temperature behavior of Oa, Jab> an< l C for 
<5.jt~0.04 and |J S |=0.09. The opening of the band gap by <5jt 
leads to the second-order transition at Tn. The inset shows 
Tn versus Sjt- See Fig. |5]for other notations. 

In summary, CO in manganites is a magnetic phe- 
nomenon. The transition from the AFM ground state 
to the paramagnetic state goes through the intermediate 
phase, in which the FM order within the zigzag chains 
coexists with the inter-chain disorder. This implies that 
the multi-critical behavior in CMR manganites [13J is the 
result of competition between different magnetic phases, 
which is only masked by spurious CO. Details of the 
phase transition at Tqo present an interesting and so 
far the least understood part of the problem. We would 
like to stress two points, (i) The transition should be dis- 
continuous, so that the change of the electronic structure 
should be sufficient to open the gap and self-stabilize the 
anisotropic magnetic state, (ii) The transition should be 
accompanied by the lattice distortion, which is a neces- 



sary precondition for insulating behavior above Tn . Both 
features are clearly seen in the experiment Q]. An inter- 
esting question, which still remains to be understood is 
how the FM zigzag chains will break. An important step 
in this direction would be the explanation of incommen- 
surate CO superstructures observed just below Too [Hi- 

I thank all member of Tokura group for the great deal 
of inspiration by numerous experimental stories about 
CO in manganites. 
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